AN ANALYTICAL INVESTIGATION OF DYNAMIC LOADS 
ON INVOLUTE SPUR GEARS 


By 

HANAMANTA G. HUNASHIKATTI 


M £. 
IH*73 



INWAN INSTITUTE OF TECHNOLOGY KANPUR 



AN ANALYTICAL INVESTIGATION OF DYNAMIC LOADS 
ON INVOLUTE SPUR GEARS 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 



By 

HANAMANTA G. HUNASHIKATTI 


to the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

JULY, 1973 



; CCNTRAL UBRAI?¥ 

, >v>, ^256E.SL 


rt? SEP '^9/3 



rnE- r‘=):> s ~ P^UN-- 



TO Mr PilEENTS 






Certified that tMs work cn ikoalytical Ijovestii^-o-tioD. 
of the i-^ynajaic Loads on Involute ^]yzT Ceo.rs” by Haiiuoanta, G* 
i-jiunD-sIiikatti has been carried out under 105^ sunorvision and that 
this lias not been siibnitted elsewhere for 51. degree* 



iissista.nt Professor 

^'ept. of Mechanical Sngi nee ring 

Indian Institute of Technology 

.Sannur 




I would like to express nj deep appreciation and 
pratitude to d'r, J. CZiakraborty for nis able .puidance and constront 
encourageraent throiagliout t'ais investigation, 

ify thanks are due to Hr. 3.3. Bhadoria, Hr* S.3. 3’nande 
and Mr. it.E. ISiare for the useful lively discussions, we used 
to have, on the subject. 

}fy thanks are also due to i:y friends, who directl3r or 
indirectly helped me in this work. 

Lastly , I thank Mr. Ifaa d-aiaan Tandey ■ f or typing the 
manuscript neatly and patiently. 


■■idlMiKIWA. 3. HCl'KiSHiraiLI'I 



CCi-TTEi-ITS 


ar'G 


CSilPTER I 

; ii^irrzDucTiai 

1 

1.1 

Gr G ne ral Int r oclu ct i ■ n 

1 

1.2 

detailed Zisciission sn.Zyncjnic Load 

3 

1.3 

Lite r a,tn re Su rvey 

7 

GEII-'TSE II 

s LEUduYSIS MD CCICUTZTICM TE'SFI?;!!! 

11 

2.1 

Elastic Deflection 

11 

2.2 

Contact iiatio 

10 

2.3 

Hesh-stiffness and Static Zef lection 

23 

2 .4 

'^nanic Response 

24 

2.5 

C onpnt at ioa-Zlgoritan 

27 

c::-L'ii-TSE III 

: LESILT'S IdZ riGCUSSICI'lS 

35 

3.1 

Static Zeflectirn 

35 

3.2 

Mes’a-stiffaess and Illative Displacement 

35 

3.3 

Dynanic Loads 

37 

CZEPTSE lY 

: GCiKiiDSIDKS ^ 

58 


lEPlKDlX A 

rJilFSEEHeES 


CCMi-ETEIi LISTItKJ 



LIST CF FIGUI33 


Pa,.rrQ , 


Fig. 2.1 

Gonfomal I-Iappinc 

ID 

Fig. 2.2 

Teeth profile - i^ctual and Mapped 

19 

Fig, 2.3 

3ertz d^efomation 

20 

Fig. 2.4 

iddvance due to relative Deflection 

20 

Fig. 2.5 

Contact -advance 

21 

Fig. 3.1 

Single Gear Tooth deflection 

41 

Fig. 3.2a 

-^ef lection-enrves for 1:1 Gear .-atio 

42 

Fig. 3.213 

d-ef lection curves for 1:1 gear ratio {■^^--dlxis- 



divisions in the path of contact 

43 

Fig. 3. 3 

Conparison nth Fnerpy Method 

44 

Fig. 3. 4 

deflection for different Gear ---otios ( 1 :1.5, 1:2, 



1:2*5) 

45 

Fi,g.3.5 

Coobined Mesh-stiffness . curve for 30P-3GG pair 

46 

Fig.S.a 

Static relative displaceoient and dynanic -response curve 

47 

Fig. 3. 7 

dynaxaic Load on Individiial Pair in the Interval 



of one Base Pitch (20P-33G .Pair) 

40 

Fig. 3. 8 

dynacic Load on a Tooth as it posses through 



the ; Snpapenent ^-vanpe . 

49 


dynanic Load' Fact or: for different speeds xfox" Gear, 



Sets with' d:i;. Gear 'Batios*/' , . ^ 

50 

Fig.3 . 1.3 

Con tact Bati o-i'n crease; '■with; ■ Load ' for,' 20P-40G' ' sets: , ■ 

51 


Fig^A^l Tooth thickness^ ' ■ , 61 

Fig,^iw2;' . ■ vGo^orcii,n'at,e ■ ^Transf o'mati OU: 


,63 



LIST T/J3LES 


Table 

1 

Happin.^ Coefficients 

34 

Tabl e 

2 

-■^ynainic Lead factors 

43 


Table 3-1136 flection Data 


5S to 



SYHCP3I3 


In tMc TJcrl:, the ieflections of teeth of gears in 
engageneiit are i.-r/estigated for different ccoibi nations of gears, 
fhe nesh-stiffness and static reiative displacenents are evalnated 
iron the deflection data. An equation of notion is fcrrjulated for 
the gears in engagement, The equation is solved by the fourth Order 
-'Unge Eutta ilethod. The dynamic loads are eTOlnated and the manimin:: 
locdsfor different speeds a.nd different combinations of teeth 
evaluated. 


are 



CSAPTES I 


niTBOBUCTION 

la INTaODUGTIGH 

Gears are essential elements of all machines. The 
efficient design and operation of most machines and instnoments 
greatly depend upon the design of gear-transmission* Gear drives 
have a considerable influence on the i^reight, sise, reliability 
and initial and running costs of the macliinos. Tb.e durability 
of the gear-drives depends upon the accurate prediction of the 
load canying capacity of the gear-tooth. Some cf the limiting 
design factors in specifying the load-capacity are : 

1* Strength : This is the load the tooth can carry 
without permanent deformation or fracture. Strength 
failure normally occurs ^tith a fatigue fracture 
initiated at or near the root fillet. 

2. Profile-durability : This is the load that can be 
carried without damaging the profile. Nomally, the 
limiting factor is the. pitting of the profile at or 
or below the pitch point. Occasionally, when 'high 
■ overloads are experienced,' plastic' profile.' deforma- „ 
tion is encountered. 

Scoring Loads exceeding 'this'' Vrating' bre'ak/dom 
:' the.' oil film'..and allow, the. 'eontact in g. surfaces to', ' 
weld together, ■; ab,rade,^',,and^,otherwis,e .^de'stroy ..t^^ 

. selves. This rating is particularly to be considered 



for highly loaded case-carburized and casehardened 
gears cr gears operating at high speeds, high tempe- 
ratures, and with thin oils. 

4. Noise : It is the result of high speeds or heayy 
loads, which excite vibrations. The vibrations are 
also due to the errors in engagement between the 
meshing teeth. 

All those factors are fui:tctions of the tooth loads and 
an accurate determination of these loa,ds becomes very important. 
These loads can be classified as : 

1. Transmitted load (tangential load) 

2, Bynanic load. 

The transmitted load is the useful force transraitted 
from one gear to another during action. It is the function of the 
external loading conditions. The dynamic load is the maximm 
inst-antaneous force acting between gears during cacti on. , The 
dynamic loads are mainly due to the imperfections in the produc- 
tion nethods ,, the material properties and imperfect positioning. 
These factors produce. a, critical error in action which in the 
final analysis/may.' produce gear tooth' overloads in /excess , .of . 
safe/ loads. The' clynanic loads seriously .affect, the / load-tra 
eitting eapacity ; of gears .' and,,, a^re ■'/also, a;' 'source , ,f or ' noise',; ,' ' The; , , ' 
data ,„obtaiiie.d, :So ,'far on the .'basis' of' theG'Tetic.al"',.'arid'' experime'ntal' 
i'nvestigati'bhs. ,'of ' dynarai'c, loads .not ,:s.uf.fi,ciont'ly..'',.'r^^ . 

"for "use ',to,,.''be,' nadeiof ^then'in,’ calculations for',,det.e;.minih^^^^^ ',/".'. 
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optinun dinensions and weight of gear-drives , Therefore 

m 

work on the dynatic loads between gear teeth is continuously- 
extending. 

1.2 SBTAIIEI) DISCUSSIOKTS ON DYNililC-LOiDS 

Gears are generally designed to mn at -unif ora speeds. 

Any disturbance vfhich will affect the unifomity induces vibra- 
tions and hence gives rise to the dynaaic loads. These distur- 
bances nay he due to the errors in geoaetrical size, fom and 
position of its several features and also due to the elasticity 
of teeth. In addition to this the effect of friction and danping 
should also he considered for accurate evaluation of dynanic loads 
The errors causing the disturbance can be classified as follows. 

1.2-1 GEOMETRICAL EEROSS : 

Geometrical errors in gears i.e. the deviation of the 
tooth profile from the desired one, the relative positioning of 
the teeth on the blank and the dimensions of teeth are due to the 
cumulative effects of the following imperfections : 

1. Imperfections in cutter 

2. Inhonogeneity of the material of the gear-blank 

3. Error in indexing -while cutting or generating 

4. Inaccurate clamping of the blank, cutter, etc. 

For study and measu3rement purposes, the errors in 

finished gear sets may be classified as ; 
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(a) Prof ile— errors : These indicate the departure of the 
actual from the design profile. The principal 
sources of tooth profile error are - 

a) Inaccuracies of the generating cutter tooth- 
profile 

h) Errors in setting of the cutter 
c) Departures fron unifcmity in the relative 
notions required for generation. 

The nature and nagnitude of these errors 
obviously depend upon the type of generating process 
eaployed, the condition of the nachine and care taken 
in setting up. However in the absence of the know- 
ledge of the actual distribution of these errors, 
for analysis purpose they are generally assuned to 
be a periodic function. This assunption is quite 
logical because any type of error-distribution can 
conveniently be represented by fpurier series. 
Especially while cutting gears, the sources of error 
for generating one tooth will be the sane for genera- 
ting the other teeth also. Hence considering gear 
as a whole, the inperfect ionswi 11 be periodic, 

(b) Pitch error ; Spacing errors are largely a direct 
reflection of the inaccuracies in the indexing 
nechanisn of the gear-cutting nachine. Pitch 
errors nay occur in isolated fom or they nay follow 
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some sort of distribution on the consecutive teeth. 

Pitch errors increase the error-in-action. 

They are a potential source of dynamic load and noise. 
The fretjuency and the rate of change of the pitch 
errors from tooth to tooth greatly affect the dynamic 
response of gear-drives. 

1.2-2 EHROa DEE TO ELASTICITY OF TUB MATSKIAl : 

Because of the elasticity of the gear materials, the 
mating teeth are deflected in action. These errors are super- 
imposed on the geometrical errors and change the error in action. 

The deflections produced under the operating loads are not 
constant during the engagement of a pair of gear teeth mainly 
because of the change in the load position and the varying flexi- 
bility of the gear tooth along the profile. Moreover when the 
contact ratio is more than unity, more than one pair of teeth 
comes into action in the meshing- range. These pairs act as 
elastic springs in parallel and the resultant stiffness differs 
from the single-pair stiffness. The result will be a periodic 
variation in the combined mesh stiffness. In addition the deflec- 
tions under load of the tooth pairs cause an increase in the 
contact ratio, the increase being proportional to the amount of 
deflection and hence the applied load. This periodic variation 
in the laesh stiffness acts as a source for vibration and hence 


the dynamic load 
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In addition to these main sources, there are external 
factors which affect the dynamic load : 

1. Torsional deflection and bean bending of gear-shaft 

2. Mounting arrangement and alignment 

3. The imbalance of rotating parts 

4. Temperature differential throiighout gear box 

5. Varying output and input characteristics. 

Actually dynamic load should be detemined for each 

and every new gear application. This way the design engineer 
will have a hack log of infomation with vfhich to compare a given 
dynamic-load calculation. If a deligent procedure is set up 
whereby the apparatus is calculated in the design stage and 
observed in service in the field and maintenance and failure 
records are actually analysed, the dynamic load effect will soon 
have as much meaning for the average gear-designer as does the 
transmitted load. 

There are, of course, many well established applications 
■where it is not necessary that dynamic load calculations be made 
since past experience has already established the amount of 
dynamic load present. There are however, several types of gear 
arrangements where it is quite important that veiy careful consi- 
deration be given to the dynamic load effects during the design 
stage. In these arrangements the omission of a dynamic^load 
evaluation would be temed poor design practicei^lfes} , ^ ^ 
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1.3 LITEiRATTJRE SUHVEY 

The load carrying capacity of gears has been a subject 
of interest to uiaiiy an investigator. The wrks of Marx and 
Guttler and Boss , Buckingham [s] , and Walker are 
the pioneering contributions in this field. They have developed 
emperical expressions for dynamic factor as a fxinction of peri- 
pherial speed. An extensive study of the dynamic loads on gears 
was conducted by the ASMS Special Eesearch Committee on geai-s 
headed by Earl Buckingham £31j. The findings of this study had 
been the basis for the later-day design practice. Similar study 
was conducted by Tuplin and Reswick . The main features 
of these studies are that the gears are replaced by simpler repre- 
sentative models, simulating the impact of the engaging teeth and 
the transient effect of the isolated tooth-errors was considered. 

No consideration was given to the actual error distribution and 
the stiffness variations in these investigations. However the 
emperical relations for dynamic factors given by these investigators 
are still useful for predicting the dynanio loads at prelininaiy 
design stages. 

Strach^sJ is the first to study the periodic excitation 
due to the error distribution and stiffness variation, which 
induce vibration and cause large dynamic loads. He considered 
unmodified teeth and formulated an equation of motion with a 
consteait stiffness for a pair of teeth in contact and twice its 
value for double pair contact. Only the forced notion was 
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considered in the solution. The possibility of loss of contact 
was not taken into account. Seman||9) included the effects of 
idealised periodic profile error but neglected the variation in 
stiffness. Eiohardson£jo] studied static deflections of unmodified 
tooth profile. Ee also solved an equation similar to Strauch with 
an addition of the clamping and backlash. Ha.rris studied the 

transmission characteristics of the gears with modified profiles. 

He drew the static transmission error curves using the tooth- 
deformation provided by Weber . These curves were used to 
determine the distribution of forces between the tooth pairs. 
Simplified error curves including the salient features of the 
actual curves were used in the analytical solution. He showed 
that the amount of damping present in the system has an important 
influence on the nature and properties of the oscillations set up. 
Resonance frequencies were predicted. R.W. Gregory ^21^ in 
association with Harris and Munro continued the experiments of 
Harris and confirmed most of the predictions of Harris. Attia|l6^ 
for the first tine measured the actual tooth deflections in action 
hut most of his experiments were conducted at speeds well below 
the resonance speed, Utagawa and Harada |,18| and Neinan and 
Settig |^13£ studied experiment ally the effects of intentionally 
produced pitch errors, Utagawa and Harada also obtained graphi- 
cal solutions without considering the damping effect, A. Hichmari4_27^ 
has obtained s Glut i ore considering the damping effect and the tooth 
■ nodif icati'on8land;Qeasured:';prbfi-le..:error, 
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In the present yrork, the tooth defoitiations using the 
ccaplex variable technique as applied to t-wo dinensional theory 
of elasticity are studied. Ustinenko fsoj, for studying bending 
stress in gear-tooth, transfomed the convex boundary line of the 
tooth profile into a straight boundary line by using conformal 
napping. G.N. Baronet |25l used a similar transform functions 
with lesser nunber of terns and using the stress fimctions given 
by Aida and Terauchi^l^ end studied the deflections of a single 
tooth. The work has been extended here for a pair of teeth in 
contact. The combined deflections including Hertzian contact 
deformation along the path of contact are studied. These deflec- 
tion data are utilized in finding the static deflection charac- 
teristics and mesh stiffness for a pair of gear in engagement. 

The dynamic equation of notion considering damping, the variable 
mesh stiffness and error in action between contacting teeth is 
foimulated. The equation is solved numerically on a digital 
computer by using Fourth cider Bunga Sutta method. 

The following nay be considered the contributions of the 
present dissertation. 

: A comprehensive analytical study has been conducted 
to evaluate the combined mesh stiffness of the gear 
p^^^ at different speeds and different conhinations of 
the teefh in the pinion and the gear. Previous 
investigations were done on the basis of either^^^^^^^^^^^^^^^^^ 
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experinentally got or assimed confined inesii stiffness* 
: A theoretical study of the increase in the contact 
ratio of the gears in action has been done. Hitherto 
only one investigator Nienann studied this aspect 
e xp e ri n e n t a 1 ly , 



CHAPTEE 2 


ANILTSIS AJSD COMPUriTIGN-TSCSNIQlDE 


2.1 DEFLECTION OF G'EJ& TOOTH 
2,1-1 DEFLECTION DUS TO BSilif ACTION : 

The tooth load across the tooth fa.ce of spur gear does not 
vaiy much. Neglecting the end effects, the load can be taken to 
be uniformly distributed across the face and the tooth-deflections 


can be analysed by using the two dimensional theoiy of elasticity. 
In the foregoing analysis, the complex variable method as applied 
to the plane theoiy of elasticity is used. 

In the two dimensional theory of elasticity, the equili- 
brium in the absence of body forces, can be represented by the 
following equations 29 : 




In terns of the airy-stress 


( 2 . 1 ) 


function 'll’ , the stresses 


can be expressed -’'.s 
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so that the ecjuations of equilibriiin can be expressed as, 

TJ) = 0 (2.5) 

Tliiis the stress function U is a biliamonic function* 

Using, stress-strain relations 




n: = 


t 




X e -f 2 

X© + 2 />• 

/ ^ V 


0 ti 
'd X 

0 V 

Wj 

d ti 

Oy 


) 


( 2 . 6 ) 


where 


e = 




£3" +<r 

X y 
E 

2 (l+ U ) • 
E 


(l + U ) (l - 2 ^ ) 

= Posion's ratio. 

Substituting for CT , and 6 in terns of stress-function 

X y 3iy 

and solving for u and V ^ to' get 


a 




a 2^ 


X 


^y^ 2 ( 


^ U 






(2.7) 


dT (f)x^ 2( X + /t ) 

Now we define a hamonic fanction Pj as ' 

’ V P = A U 


a* u 


(2.8) 
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and Q-as conjugate of P, so as to satisfy the Caucliy-Eeinann- 

conditions 


d? p _ 0 Q 

0 X ~ ~~'Qj 

B p _ B q 
c^y 0 X 

jSron these functions vre define 
f(E) = P + i 0 


(2.9) 

conplex function f(z) such that 

( 2 . 10 ) 


Tfliere z = x + iy 


This function is liclonorphic in the region 


get 


Defining a new ccopiex function, 
9^(z) = P + iq = I ^ f(z) dz 



Q = 


- 4 



B q 


of the body. 


( 2 . 11 ) 


( 2 . 12 ) 


Then equations (2,7) can be expressed as, 



2/U 


4t 

0X 


_ ^ ^ TJ ^ 2(X -^ 270) ^ p 

B ( >. + ^) 

_£Le , 2. (A* zH) Oo 

Qy^ (X+/W.) S)y 


(2.13) 


Integrating and putting the constants of integration equal to 
zero, as they represent the rigid' body notion, we get • 
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2 yU- u = 

2 T = 


Qx (V ^ ) ^ 

j0 U 2( X+ 2 M) 


Not/ froa 'tLo •^ofini'tioEi of It caa bo sIict/ii t.Iia'fc, 

'(U _ ps - qy) = 0 


i. . U = + qy + 


In the conplex fom, 

U = Se£a y- (a) + ^ ( 3 )j 

xThere (z) ' 

By using the , conjugate property 

20 . S (z) + z + Z (z) 

Taking the pa,rtial derivatives and adding 
C) TT g 


u . 3 
^ " -37 


0+2 ^ -( 2 ) + ^ (z) 


liiere 




dX 


(2.15) 


(2,16) 


(2.17) 


/o -f O ^ 

\cu.ia} 


(2.19) 


5hen the displciceiients o.ro by 

,2 ^{u + 1 V ) . (- Xfp . X^)- 

= (3 - 4 l?) ^{ z)-z ^{ z )~ 

In the body .BpacGv'iSdtransf orn siripler , .region by ', 
t'ransfo'matiGp'.fnnetlon 


z = 


In the transfomed co-ordinates, the displacement z-plane, we 
can write as 

2/U.(u* i v) . (3-4!/ ) ^ ^(S) -Jf(5 ) (2.2! 

The resultant defomation is 
/^2 2~ 

+ V (2.22) 

2.1-2 TEANSFOmiATION FUNCTION ; 

The seoi-infinite space 's‘ in z-plane having a boundary 
of tooth shape is transformed into the infinite space 'B* having 
strax^t boundary line in z-plane. To achieve this Ustinenko 20 
and Aida 19 have assimed a transfomation function similar to 
the one given below ; 

..4 . _ :■ 

^ (2.23) 

n=l X + b 

n 

This function satisfies the two necessary conditions of confomal 
napping, namely 

, V: 

: ■ : : ^ (.5 :: ^o: : ; - ^ 

where ir 

Here u and v do not represent the displacements but they represent 
the coordinates of the mapped space as shown in Fig. 2.1. 


16 


2.1.3 STHESS FDNGTIONS : 

When a concentrated load P is acting on the point (0,C5<) 
on the boundaiy of a semi-infinite plate (Fig. 2-l) , and foms 
an angle with the x— axis in z— plane then the airy— stress function 
U is expressed hy 19 J 


^ “ 2~^ ^ log(z-i X )e“^ ^ -(z + i -X )log(z-i(5t )e^(2.J 

Then from (2,18) 


^ 2 TT ^ 0<.)log (z-i o( )+ i pC e^P log(z-i 

t (2:2 

If7(z)=^^(z) 

the n we can write 


‘^1 (=) - 


2 /r i-" " y+e » - e 7- / (2;2 

Let ^ be a point on the boundaiy, 's' be the length along the boundaiy 
curve, the boundaiy condition to be satisfied is 

¥'(-r ^ • - 


()*Co{ f ) ^ (-f )- 1 / [s^(sKWy(s)] . ds (2. 


28 


To satisfy the condition the stress functions are taken in the fom 


- 


such that are having singular point on the boundaiy 
and ^ are having no singular point and — 6^ 0 as 


i2l2 
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Then from (2,26), and (2.27) we get 


log ioL) 


{^- i(A) 


^ (-i O^' ) 

6t)^(cX)C^ + /I-) 


(2,30) 


®y ^Pplyi^ig tile boundaiy condition and integrating over the 
houndaiy, 

^ 5^ fo ^ == 'Ti L £^, (("-JJ ) "f . 


-f (5) - 103 - io<) - " r 

' (2.32) 

Fron the nature of the transformation function CO ( ), the 

poles stand at points ( — h^) and hence using the residue of each 


27ri 


't^(_ f)f( 


pole, we have. 


1 “ ("n) 

, / , a / (h ) 

2,4 DEFOBhlTION BUS TO HEBTai/J? CCNTiiCT LOAD 


(2.33) 


(2,34) 


7fhen a pair of teeth cones in contact and cariy load the 
contact profile suffers from Hertzian defcmation at the contact 
area. The half width of the strip of contact between the meshing 


teeth is given by Hertz as, 


^ = (^) ( 


TT ' ' (l-f f2 


iC ,1^-^ % ) 


2D 1/2 


(2.35) 



la 


where f ^ fire the radii of curvature of the surfaces in contact 


1’ 1 2 
and W load in lb. 


The resulting deforoation of the neshing teeth is given 


by Vfeber£5j as 

2 W (l- j 2 h 


"H TT 


log 


1^1 


L 


b 

2 h, 


2 (1- U ) 


2 W (1- Up 

.J. 


J 


E. 


log 


2 


2{i- i;„) J 


] 


(2.36) 


The addition of individual deflections of each tooth and 
the combined contact defo mat ion yields the total deformation fron 
which the total stiffness of the individual pair is detemined.. 


= Zg + Zp + Zjj (2.37) 

Zg = Gear tooth deflection 

Zp = Pinion tooth deflection 

z„ = Hertz deflection 
H 

Zp = Total deflection 
2.2 . CONTAGT-EfUIO 

For a gear-set with conta.ct ratio noro then unity, when a 
pair cone in contact there is a certcin amount of nis ^aatch 
because the pair in engageaent has suffered a defGmation’e’(Fig,2— 4) 
TlK inconing profile ’zd lags behind by the sane amount ’e’ and-', 
hence the contact starts at a position ' s' earlier than the 
theoretical position 's'. Similarly the disengagement is post- 
poned, Effectively the contact ratio is thus increased. This 
has marked effect on the conbined mesh stiffness. The exact 
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positions of engagenent and disengagenent can be calculated 
by the geonetiy as foiloTTs 25 • 

Fron triangle BS4 in Fig, 2.5 
Lma . Sin-1 ( gl Sin (90 ^ 


10 


OC = 180 - (90 + J3^) - £bEA 
Fron triangle AGE, 


CE = (^0 ^ cs . 0 Cos « )i/s 


>1 -X* I 


... JL 0. 

^0 


ihen 


DC . (C^ 4. - 2 C Goscx ) 

, . c.- -1 Po 

fA^Sxn (- 3 ^) 


i/2 


^ =Ai-°g®g-/3 

But by involute geonetry 

^ « inv - inv 

Tifcere Cos"^ (Eg^/^C) 

and inv “ 'fcan A - X . 

Tbe symbols are to be taken to repiesent tbe quartities as 
abown in:Fig,v2.5. \ 


■ 


(2.38) 

(2.3f) 

/ 
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3!Si.e adramee angle determining the first contact position is 
obtained by the solution of simultaneous eq_uations (2^38) and 
( 2 . 39 ). 

4 • 3 I^iBSIi STIFFNESS m) STATIC DIS?IACEI#]MT 

Two successive .pairs of 'seeth come in engagement with 
a phase difference of the base pitch. The displacement charac-- 
tori sties are periodic with a period of the base pitch. The 
combined mesh stiffness is evaluated by placing the individual 
stiffness curves at a phase differenceof base pitch. The actual 
engagement and disengagement point for any tangential load are foioid 
out by using the analysis of the section (2.2. ) The stiffness of the 
individual pair in the region beyond the theoretical region is 
obtained by extrapolation* 

From the combined mesh stiffness and the profile error, 
the static displacement is computed as follows: 

W 4- ei kj[ 4- ©2 k2 

xs = — 

■where 

Yl = tangential load 

e^^ = combined profile error for the first pair 
eg = corabined profile error for the second pair 

- instantaneoua mesh stiffness of the first pair 
kg = instantaneoiis mesh stiffness of the second pair 
k =instantaneons combined mesh stiffness 
Profile error here means the combined error due to geometrical 
error and total elastic deformation ’e'. 
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2.4 DYMMIG RESPONSE 

In the dynamic nodel a pair of gear in engagement is 
replaced by two masses connected by a spring with time dependent 
stiffness; the masses being the equivalent masses of the gears 
placed at the base circle and the stiffness is equal to the combined 
mesh stiffness. The equation of motion for the system can be 
formulated as follows : 

Symbols : 

- total load 

P - static load 

s 

T - shaft torque 

~ frictional torque, . 

I -- mass moment of inertia ' 

M - mass at base radius 

K - stiffness 

Por the dynamic equilihriim of individual gear, 

■', " ■ ■ T' • . 

Similarly, ©g = (P^ " 

■ "^2 2 ; 

and ©2 are having different signs, 

But the combined displacenent is given by 

= - hi “2> 


X 
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Therefore 


« • 

X 


- (r ftg) 


Substituting for and 0^, we get 


(2.43) 


X - 




f2 


(?m - P ) { ^ + — ) - :: + 

12 12 


) 


T 

i*C* X 4- ( — ii ^ 

h 


But 


1 


f2 




■) + P. (~ + ~) = P + 7^) 


(2.44) 


Xq 

-2 


so putting 


M 


m/ 


and since the second term represents the frictional resistance, 
replacing it by the damping term f i we get 


*:^ + f i + 


11 . - ^ 

M ' ^ M 
e e 


(2.45) 


f is the damping coefficient. 

If e^ and are profile errors, then the total force 
is given by 


Pj, = kjL (x - e^) + kg (x - Cg) 


(2.46) 
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Then equation of notion hecones 


k, 


X + f X + ^ (x - ep + 2 


^o(x - e„) 


¥ 


M 


Sinplifying 


.. . . ■*' -IJ ^ ^ 

X + f X + — X = — “ 


U 


M 




By putting = :^ = 


M 


Mg 

^2 

M 


^2 


xs = 


Pg ■*• \ ®2 

+ kg 


xs represents the static displacenent , 
Then we get 

2 


X + 2tU^ X + of X = 6u 




(2.47) 


(2.48) 


where = ratio of actual daoping to initial damping 

= f/2 / k Mg 

The load carried by each pair is eYaluatecI by utilizing the 
principle of two -vSprings in parall^ due care being' given to 
the, ; errors " The to qui red .equations ca-n be shoim to'/be 


P. = k. (x - e ) 


k2 (x - eg/ 


(2.-49) 

(2.50) 
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2 . 5 COMPOTATION-ALGORITH 

' ■"'■■■ ; 

Once the analytical formulation to fevaluate the sta,tic 
deflection, mesh-stiffness, the error in action, the relative 
displacement of the pinion and the gear, the dynamic response 
and dynamic load of the gear-pinion system is complete, these 
characteristics are computed for different sets of the gear and 
the pinion. The whole procedure is explained in steps as follows 


Static-deflection-calculation : 

Step 1 The mapping coefficients for the pinion and gear are 

calculated from the values of the coefficients given by 
Borenet£ J. Any intermediate value is found hy inter- 
polation, Using these coefficients, the transformation 
function is evaluated and the corresponding coordinates 
of the profile and central-line of the gear rooth in 
z— plane and -plane are calculated. The divided 
difference tables for the respective co-ordinates are 
evaluated. ■ 


Step 2 Frcaa the gear geometry the relative positions of the 
centres of the fixed system and tooth ^stem of 
co-ordinabes are fixed. 

Step 3 The engagement range for the set under censideratiGn 
evaluated. The path is divided into 100 equal 
parts and the radii of the curvature of the profiles 
at contact at these points are calculated. 
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Step 4 -it each of the successive points considered above, the 
co-ordinates of the contact point and the central 
point, at which the path of contact cuts the central 
line of the tooth are calculated. The corresponding 
tooth system coordinates are calculated by coordinate 
transf ciriation . Then the coordinates of these points 
in the napped plane are foimd by interpolation. With 
the load acting along the contact line at the contact 
point the deflection of the central point in the 
z-plane is calculated by making use of the equation 
( 2 . 21 ). The stress functions contain their derivatives 
at pole points(bj^). These are evaluated by differen- 
tiating the stress function and substituting the values 
of bj^. The resulting four equations along with their 
conjugates form the eight simultaneous equations, which 
are solved by using the Gause-Seidel elimination tech- 
nique. The bean deflections of the wheel-tooth and 

pinion tooth are calculated separately. Inowing the 

radii of curvature of the profiles in contact, the 

Hertzian deformation is evaluated by using the Webber^ s 

nodified f omula ( 2-36 ) . The sun of all these compo- 
nent s Ogives: the total deflection, from which the ; > 

stiffness of the tooth pair in contact at the position 


can be calculated. 
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Il'«3h-3iif fnoss i>ynr.’iic Itesponse 

Step 1 F rcn the lel'lccticn data the stiffness characteristics 
of a pair cf teeth, fron the beginning of the contact 
tc the and of contact for the gear set chosen are 
evaluated . The stiffness curves for the consecutive 
pairs arc place 1 at a phase difference of base pitch and 
the ccnbinecl nesla stiffness data are calculated. The 
increase in contact ratio is taken care of by finding 
the actual beginning of contact and the end of contact. 
The individual stiffnesses in this region are calculated 
by extrapolation. In calculating the actual point of 
contact, the tooth nodifications and nanufacturing 
errors are super-iaposed on the error due to the 
deflection of the pair already in engagenent. 

Step 2 The error in action ie n&de of tt;o conponents namely 
the difference in the deflections of the pairs in 
contact and nanufacturing errors. In the conjugate 
contact- region, the deflections of the two pairs are 
ooual and hence the errors are nade of only the nanufac- 
turing errors. In the non-conjugate action-region, the 
error due to the differential displaceaents are 
evaluated by the georaetry f or tip contact. 

Step 3 ^ Thfi: deflectionr-'Curves .are -evaluated for: t static 

: loading conditions 
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Step 4 The equation of notion (2,49) is solved by the classical 
^inge— Kutta Method by using the stiffness and static 
deflection data at regular intervals. The total nunber 
cf intervals taken arc hundred. The corresponding tine 
interval is calculated by finding the tine required for 
covering one pitch, depending upon the speed. 

Step 5 The dynanic loads carried by the pairs are calculated 
by using the equations (2.50), The naxinun values for 
each pair and the corresponding position are found. The 
greater of the two naxinun values is the greatest dynanie 
load and the dynanie factor is calculated. 

The dynanie factors for different conbinftions of gears 
and speeds are evaluated. 



COMPUT.Aa'ICWAT. 

i'UiGOEITHM 
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DYMMIIC BE3P0HSS ^ 

T ^ 

INPUT 

Gear details 

Individual stiffness data 
Profile error, Pitch error 
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— 


Call %-nanic Load at each 
position. Find individual 

uoxiniin and owr iifixinmn 

%naaic Factor 


C§.) 


iBDITIOi’LiL 
BOLTINE POil PITCH 
EffiiCi: 



Error in action., Mesh-stiffness, 
Displaceaent for Two Mesh periods] 


Solve the Synanic equation for 
Two Mesh-periods 
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TISIS 1 


IIAPPING GOEFFIGIEOTS 
Pressure ilngle 20® 


Teeth 

Coeff . 

20 

30 

40 

60 

80 

150 

b, 

.005 

.005 

.006 

.005 

,005 

.006 

bg 

.025 

.025 

,025 

.023 

.023 

^022 

^3 

.127 

.127 

.127 

.125 

,124 

.124 

^4 

.490 

.487 

.483 

.470 

,455 

.450 


.0024 

.0024 

.0024 • 

.0024 

.0024 

. 0024 

ag 

.017 

. ai 6 

.015 

.013 

.012 

.012 

^3 

.130 

,127 

.124 

.115 

.112 

.110 


.220 

.250 

.275 

.325 

.335 

.340 
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3. ISS11.TS DISCUSSIONS 

The results of conputation gind tixe discussion of the 
results are given below. Sene typical gear tooth conbinations 
are taken. The dianoti-al pitch of the teeth is 4/inch in all the 
cases. The naterial cf both the gears is steel having S to he 
equal to 30.10 Ib/in*' and to be 0,3. The teeth have standard 
involute fora and the pressure angle is taken to be 20 degrees. 

3.1 COIBINSD STATIC DEFLECTIONS GF THE CONTACTING PAIR OF TEETH : 

The static deflections of a pair of teeth in contact are 
found by using the complex variable nethod equations 2,21-2.34 . 

The following gear tooth conbinations are taken : 20P-30G, 20P-40G, 
20P-60G, 30P-3CG, 40P-40G, and 50P-50G, 60P-60G and 70P-70G. 

XP-YG neans that the pinion has X teeth and thegear y teeth. Sone 
of the results are shown in Figs. 3.1, 3.2, 3.3. There is a 
gradual reduction in c(xibined deflection as the contact point 
approaches the pitch point and'Tt 'iaoreaseE as the point noves 
away. The deflection characteristic is sinxlar to thooe we got by 
energy nethoil . The exact fom of the conbined deflection 

changes with nu-r.ber of teeth (reduces with increased nunber) and 

nomal load on gear tooth. 

3,2 MESH-STIPFNISS AND EEhATIVE STATIC DISPLACIMINT OP THE TEETH : 

For a particular pair of teeth (one of the pinion and the 

other of the .gear)/iu,n;,ataet:,:;:;it/i8:;^ : 

characteristics will be droping down from approxioafcely the pitch 
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point. As t!ic' c iitnct ratio of the spur gears, in general, is 
quite appreciable {ab.ut 1.7), for cost of the contact path two 
pairs will h*; in contact* Tlie resultant coribined neah— stiffness 
of the fjervr systf';^ fr.r different conhinations of teeth is shown 
xn 3*ti, Xt is soeti that as in the neighbourhood of the pitch 

point only one pair r=f teeth is trrasnitting power and because 
of the cosnbinol effect, the stiffness is least in this zone. At 
the transition peried when the contact is spreading from two pairs 
of teeth t j one pair or when the contact is shifting fron one 
pair to two pairs, there should be a sudden change in the conbined 
□esh-stiffness. But this effect at the transition stages sonewhat 
snoothena cut because of the change of the contact ratio due to 
deflections of the teeth. In this dissertation the effect of 
inpact has not been considered. 

As the variation of the tooth deflection for different 
angles, of rotation of the pinion is not constant, the transaission 
error also varies with the rotation of the pinion. Transaission 
error can be cxperinentally neasurod by narking the gear and 
pinion under no-ioad position at a coinc-idental point and observing 

the relative displaceaent 3 of these a arks when the gears are in 

notion under load . The transtiission error cuive for 

20P-30Cr is given in Figure il. 
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DYX^l’^IC HESPONSS 

fhc bcundaiy value problea (equations 2.49 and 2,50) 

-with danping factor of ,1 of critical, is solved by assuning the 
initial values of the displeceaent and velocity as 

V, = (Xs- - XsJ/ AT 

At the interval size. The tine T of one cycle, i.e. the tine 
needed for the contact point fron the pitch point of one tooth to 
that of the next tooth is divided into hundred equal divisions. 
Each of these divisions is equal to ^ T. 

The final values are conpared with initial values and the 
fcllowiji {5 convergence criteria arc applied. 

-iML 1 ^ .001 

"^,31 

( ^ .01 

The convergence to the desired solution has occured in about *-5 
iterations. The (fynanic response at 500 BIM for 20P-30G is shown 
in Pig. 3.6. It shows peak values whenever a pair enters or 
leaves the contact zone. In addition to these high peaks, 
mailer peaks also oociir in between because of the vanatiGn of 
individual tooth stiff ness,.^ , ; '-, 

' ; ' / '-^The- dynaiaie'^ 

la- Pig,: 3.7, ;^;It;^ahc»r 

load on a pair ^ has a ^ 
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investif'at, rs after conducting oxperinents observed that near the 
pitch pciiA the pitting phenonenon was nost proninent and it is 
accepted tiiat when the contact is near the pitch point, the gears 
exi'erience the worst loading condition. This has been substantiated 
by the presctit inirestigation also. The worst dynanic loads occur 
when the c vntnet is near the pitch point. Location of the exact 
spot where worst loadincj: occurs is not feasible either by experi- 
raental or theoretical noans. The author feels that as the kinematic, 
geometric and load characteristics do not change appreciably in 
the neighbourhood of the pitch point, the pitch point itself can 
appr,.xia{±oly be taken as the point representing the worst. 

In Fig, 3,7 dynanic load of an individual pair of teeth 
in contact ia shewn a^jainst the location of the contact point 
frcci one pitch point to the next one. The sane result is given 
in Fig, 3,8 to show how the dynamic load varies on this pair of 
tooth as it cmes in contact aind leaves contact. The dynanic 

load factor considering the naxinun peak lead has been conputed 
and is shown in Fig. 3.9, • As expected the factor varies 

with speed. It increases upto the critical speed of the set. 

Beyond this it cones down and tends to be constant as the speed 
increases. The dynanic load factor for speeds well above the 
critical speed does not vary mch for different gear sets. The 
naxinun dynanics-load factor varies for different sets. The 
hipest value ie about 

those found experinentally and predicted by sone investigators. 
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Kv.'ihlci I I'f J Cl uclui-ci fron his test jrssul'ts that irrssDcctivo 
c f err. r, the lyna^iic load never exceeded three tines the design 
I jivi, The values calcul <ied by Buckinghans equation | 11 '] are 
shawii in Table 12. They are also telov the critical value of 3. 
Harrin Ktr.tes that f *' r nolified tooth— gears running at the design 
load, the fact' r will have a nnxinim value of two. There has not 
been universally accepted rating formula for the dynanic factor. 



Tx^LE a . 




iyN/i!JIC LO/iD PACTOna 
Tangential Load 1000 lb. 


S.No. 

1 Goar Sot | 

Speed 

EtM 

{ %nanic 

1 Load (lb) 

1 %nanic Factor 

2.86 

1 

30?- 30G 

500 

2859 

2 


1000 

1200 

1.2 

3 


1500 

1559 

1.559 

4 


2000 

1666,0 

1.66 

1 

40P-. 40G 

600 

1665 

1.6 

2 


1000 

1229 

1,229 

3 


1500 

1350 

1.358 

4 


2000 

1650 

1.650 

1 

6QP^ 50G 

500 

1268 

1.268 

2 


1000 

14060 

1.4 

3 


1500 

1634,6 

1.63 

4r; ■, 


2000 

1639 

1.64" 


DEFLECTION -IN lOOO 


0,7999 

0.6999 

0.5999 

0 . 4999 ” 

0. 4000 

0.3000 

0.2000 

O.lOOO 

■ 0.0000 









Fig. 3 3. Comparison of energy and complex variable method 
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TilBia 12 

mmMIC LOJQ P/iCTOIJS by BUCim-JQlLM 
EQUATION [3^ 


Design Load | 

w 1 

11 

Dynanic 1 

Load 

i 

Ilatio 

1 / TJ 

625 

1567 

2.52 

725 

1465 

2.02 

2080 

4580 

2.2 

2939 

8326 

2.84 

4160 

9770 

2.35 

5270 

12195 

2.31 
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4, CONCLUSIONS 

4. 1 The following are the nain ccnclusions drawn out of the presen 
imestigation ; 

a) The worst loeding position of the gears is near the 
pitch point. 

b) The naxinum dynamic load factor is 2,86, 

c) The dynamic load factor is different at different 
speeds. It attains the naxinun value near the resonance 
frequence of the systen and as the speed increases it 
decreases. At higher speeds the dynamic factor for all 
gear tooth conbinations is apx:-roxinate equal, 

4.2 SCOH: FOS FUi^THBE 1701$: ; 

The present investigation has the following features in 
c caparison with other work in this field ; 

(a) It is the first attenpt to evaluate, through theoretical 
analysis, the conbined aesh stiffness of the gear systen. In 
earlier work, the nesh stiffness characteristics were either 
chosen arbitrarily or an estinrhe was made experimental ly, 
Nakanura 23 chose a square wave pattern and Seireg 25 
analysis with a trapezoidal pattern. 

(b) In almost all ihe in-rostigations, the transnission ratio 
was taken to be unity. In this case some of the character- 
i sties are idential about the central positions. If the 
numbers of teeth of pinion and gear are different, this 
symetiy about the central positions is destroyed. In the 



present inTOstigation gear c rot inat ions with unequal 
nunber of teeth have also been analysed. 

However the present investigation is far fron the ideal one. 

S one significant assimptions are made here. In order to 
refine the results obtained, this work should be extended 
in the following vray : 

(a) The daizping factor should be detennined accurately and 
the analysis should be done with the refined value of 
the daaping factor. 

(b) S one provisions should be nade to take into account the 
frictional forces and the impact of the teeth. 

(c) Though provisions have been nade in the present analysis 
to take accounts of ‘tire" effect of the pitch error and the 
profile error, no ccoputat ion is perfomei with these 
errors. It is suggested that the dynamic factors should 
be computed in the presence of these errors. 
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p = i 


0 . 


(b) THE INCLIN.4TI0N OF THE NOmiM. TO THE LINE 

INCLIN /JION TO TH E BASE-RADIOS 

, (f )- 

= tein^ - 

TT 
TT 


Q. 




o 


p> 


(4) 

(5) 




= radius of cuivature a of the profile at the position 


(c> CO-OBDINilE TSa^FOmaTION 

The co-ordinates of the contact point in fixed systen 

(") 

The co-ordinates of the intersecting point on the median 

Q (- f - h) 

where h = (tan sj-- - tany3 ) 

The transformation to the tooth-hase-system is carried out 
according to, 


X- 


X 

s 
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= M M 

or rs ■ 

1 

Ys 

1 


1 
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or 
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rs 


Where 
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-X 
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0 
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Cos j3 
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